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STABILITY REGIMES IN A ROTATING 
QUADRUPOLE FOCUSING ACCELERATOR 
I. INTRODUCTION 

A number of recent high current accelerator configurations utilize strong focusing 
fields. These fields, consisting of a stellarator field (or rotating quadrupole field) and an 
axial guide field, increase considerably the energy mismatch tolerance of the device ancl 
provide confining forces against the beam space charge forces. 1,2 Two such devices are 
the modified betatron accelerator 3 and the spiral line induction accelerator (SLIA). 4-9 
The recent addition of strong focusing to the modified betatron at the Naval Research 
Laboratory has allowed that accelerator to successfully accelerate a 0.5 kA beam to 12 
MeV for ~ 35,000 turns in a vacuum chamber with resistive walls. 7 The SLIA, to be 
constructed by Pulse Sciences, Inc., will utilize the strong focusing for transport along the 
curved sections of the beam line between the accelerating cavities. 

The use of strong focusing fields has a potential difficulty in that they can lead to 
various types of beam instabilities. 8,9 It has been suggested, for example, that the rotating 
quadrupole field may act like a wiggler field as in the free electron laser. 8,10 Of interest in 
this study, however, is the three-wave instability in which the externally imposed rotating 
quadrupole field interacts with the transverse motion of the beam centroid to excite a 
transverse-electric (TE) waveguide mode. 9 Approximate stability conditions for such strong 
focusing fields have recently been found. 11 

The purpose of the present study is to obtain detailed analytical conditions for stabiliz¬ 
ing the three-wave instability. Here we analyze the dynamics of electromagnetic waveguide 
modes and the beam modes associated w' n a relativistic electron beam propagating under 
the influence of external strong focusing fields. The dispersion relation for such a system 
is generated in Section II below. In Section III, we show that this dispersion relation, in 
the limit of zero beam current, gives analytical conditions for the various stability regimes. 
We find a) two physically distinct three-wave unstable regimes, b) two three-wave stable 
regimes and c) a regime in which the particle orbits themselves are unstable, irrespec¬ 
tive of the electromagnetic waves. The orbit unstable regime is a recovery of an earlier 
result. 1,2 In addition, we obtain algebraic expressions for the growth rates in each of the 
three-wave unstable regimes. This is done with the simplified dispersion relation where 
the right-hand (RH) and left-hand (LH) circularly polarized waves are decoupled in the 
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low current regime. These are found to be in general agreement with the approximate 
expressions given in Ref. 12. Numerical solutions of the dispersion relation are presented 
in section IV. The analytic expressions of the growth rates and stability boundaries are in 
good agreement with the numerical results of the dispersion relation. 

II. MODEL 

In our model the external fields consist of a periodic rotating quadrupole field and a 
longitudinal magnetic field (see Fig. 1). The electron beam radius is assumed to be small 
compared to the waveguide radius. In the equilibrium position, the beam is centered along 
the axis of a circular waveguide and the beam electrons are assumed to be monoener- 
getic with zero transverse velocity. Perturbing electromagnetic waveguide fields cause the 
beam centroid to develop a transverse velocity and become displaced off the z-axis. This 
displacement of the beam centroid amounts to a transverse macroscopic current which, un¬ 
der certain conditions, further excites the electromagnetic field. The displacement of the 
beam centroid also induces image electric and magnetic fields on the wall of the waveguide. 
The beam centroid motion, in our model, is governed by the following fields: i) rotating 
quadrupole field, ii) longitudinal magnetic field, iii) electromagnetic fields, and iv) induced 
image fields. 


Wave equation 

Before developing the orbit equations for the beam centroid we first derive the wave 
equation for the electromagnetic fields. The electromagnetic fields are represented by a 
vector potential given by, 

A = A + (x,y,z)e~‘ w< + A_(*,y, *)«“*"* + c.c., (1) 


where A + and A_ are complex amplitudes associated with the right-hand (RH) and left- 
hand (LH) circularly polarized waves in a cylindrical waveguide. In Eq. (1), is the radian 
frequency and c.c. denotes the complex conjugate. 

The wave equation for A is given by 


(—2 <9 2 1 d 2 \ 4a- 

+ dz 2 ~ ?dt*) A “ ~T 3 ' 


( 2 ) 
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where Vj_ is the transverse Laplacian, and J is the macroscopic transverse current asso¬ 
ciated with the beam centroid. Substituting (1) into (2) and operating on the result with 
dtexp(iut), in order to select the correct frequency dependence, yields 


L( A + 4- A_) = -- e iwt J, (3) 

C J 0 2tt/w 

where L = V^_ -I- d 2 /dz 2 + u 2 /c 2 . To obtain the RH and LH polarized components of the 
current, we equate the e+ and e_ components of both sides of (3), and find the following 
wave equation, 

4— r2ir/uj j a 

L( A + +A.).(2e T ) = - T J — (J ex TiJc V )e iui , (4) 

where e± = (e* ± ie y )/2 and e x and e y denote unit vectors in the x- and y-directions. 
The current density associated with the beam centroid motion is 


3(x,y,z,t)=-\e\X h 

6(x - x e (z 0 ,t))6{y - y e (z 0 ,t))6(z - z e (z a ,t)\. 


f dx c {z 0 ,t) , 
dt e * 


dye (*o.O . 


(5) 


where A* = nj7rr£ is the number of electrons per unit length, is the beam density, rj 
is the beam radius, x e and y c denote the transverse coordinates of the beam centroid, 
z e denotes the axial position of a cross sectional slice of the beam and z„ is the initial 
position of the slice: z e (z 0 ,t = 0) = z 0 . In the small-signal or linear regime z c and y c are 
proportional to the electromagnetic fields and (5) can be written as 


J c (x,y,z,<) ~ - |e|A fc £(x)$(y) 

(dx e {z 0 ,t ). , dy c (z 0 ,t) . 


/x- 


at 


at 




where v„ is the axial beam velocity. Substituting (6) into (4), we obtain 


4-r r2ir/u> j. poo 

I(A+ +A_)-(2e T ) =—v6(x)6(y)j' — 

(dx e (z 0 ,t) .dy c (z 0 ,t)\ 


dz„ 

— z„ — v 0 t)e' 


( 6 ) 
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where u = (|e|/m 0 c 2 )/&//? 0 ~ It,[kA)/\7t3 a is Budker’s parameter, If, is the beam current 
and j3 0 = v 0 /c. Upon carrying out the time integration in (7) we obtain, 


L( A + + A_) - (2e T ) =^uS(x)S(y)e^^ 



d 


2 ttv 0 /u dt 


(x c Tiyc)e- lUJZ '’ /v ", 


( 8 ) 


where x e (z 0 ,z) = x c (z„,t = {z - z 0 )lv 0 ) and y c (z 0 ,z) = y c (z 0 ,t = (z - z 0 )/v 0 ). Since 

our model is spatially periodic with period u)/2ir, the limits on the integral have been 
changed to 0 to 2ir/u/. 


Beam Centroid Orbit 

The configuration of the rotating quadrupole and the axial magnetic fields are shown 
in Fig. 1. The total external magnetic field, B eit , consists of the rotating quadrupole field 
and the uniform longitudinal field and is given by B ext = (B qx ,B qy , B zo ) where 

B qx = -B q k q (x sin k q z - y cos k q z), 

B qy = B q k q (xcosk q z + ysinfc,z), (0 a - c) 

B zo = B 0 . 


In Eqs. (9a-c), B a is the axial magnetic field, B q is the magnetic field of the quadrupole, 
k q = 2?r/A, and A, is the period of the quadrupole field. The representation for the 
quadrupole field in Eqs. (9a,b) is valid near the z-axis, i.e., (x 2 + y 2 ) 1/2 << A 7 /2t. 

In addition to the electromagnetic fields given by (1) and external fields given by (9), 
we have induced fields due to the displaced beam. These fields are produced by the image 
charges and currents on the waveguide wall when the beam is displaced off the z-axis. For 
a circular, perfectly conducting waveguide the induced electric and magnetic fields near 
the z-axis are 


m 0 c u / 


(ion) 


4 


(10 b) 







where r g is the waveguide radius and we have assumed (x 2 4- yl) l/2 << r g . The motion 
of the beam’s centroid under the influence of the fields in Eqs. (1), (9) and (10) is in the 
linear approximation governed by 


7 0 m 0 c \dt dz) 

~r~ 5 “ - + v 0 k q $l q (y c cos k q z - x c sin k q z) - v 2 0 k]y c 


= (116) 

where £l 0 = \t\B 0 j^ 0 m 0 c is the relativistic cyclotron frequency associated with the axial 
field, fl q = \e\B q /j 0 m 0 c is the relativistic cyclotron frequency associated with the rotating 
quadrupole field, = (1 — &l)~ 1 * 2 is the relativistic mass factor, k s = {2z'/(/3 2 y 2 r^)) 1/2 , 
z = z 0 + v 0 t, and A x , A y are the x and y components of the vector potential given in (1). 
For the purpose here, we will approximate the exact expressions of (1) by the fields on 


— A±(x = 0,y = 0, z) = a±e tk±: e±, 
m 0 c* 

where k± are the axial wave numbers. 

The orbit equations for the beam’s centroid can be written in a more convenient form. 
Setting £ = x c + iy c , Eqs. (lla,b) become 


K ±. 

\dz 2 ° dz 


k) U - K q k q e ik '‘C = F 


where K 0 = n o /v 0 , K q = Q q /v 0 , and 




+t(fc_ — (jj/v 0 )a-e ’ 


i(fc- — uj/v 0 )z iu)Z n /v„ 1 


It is convenient at this point to introduce a transformation from the quantity £ t 
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Substituting (14) into (13) transforms the beam’s controid equation into 


where K x = A c - Jfc„ K\ = (A 0 - k q /2)k q /2 - k 2 , A 2 = K q k q and G = Fexp( -iik,r/2). 
Equation(15) can be written in the form 

+ (2 K\ + K\)~ + (Ki - K‘)j i = H, 116a) 


where 

B = [J? + iKl lh + *»] a + K ' G ' 

= *- >k ' zn + UK, - *,)£ + ((A', - k,/2)k,J2 + A'|)j F 111)41 

+ AfFV*^ 2 . 

Substituting (13) into (16b) yields 

H = jy [(#; ~ k q /2)D+a+e~ iK + z - (Kl 4 - k q /2)K 2 a'_e ~ iK -*] 

— ~z ~— [(A_ - k q /2)D-a-e lK ~ z - (K+ + k q /2)K 2 a + e lI<+z ] 117) 

PoTo 

where A ± = k± — w/v 0 ± k q /2 and D± = A'i T K X K± - A'|. The particular solution to 
(16a) is 

f = (a;a;e" iiC + I +/5:a‘_e- ,A:lj ) e - ,a,J " /, '‘’ 

+ (a_a_e ,i: ^“* -4- l3+a + e' K ^ z ) e "~ Zni '"", 

where 

a* = -j-^D±{k ± — <jj/v 0 )/R±, 

0 ± = - Kla ± /D ± , 

A± = A'£ - (2A\ 2 + A' 2 )A'i + A', 4 - A', 4 
= (K 2 ± ~ (d 2 + d 2 )) (A'I - (d 2 - d 2 )) , 
d\ = A, 2 + A'f/2, 

d 2 = ((A 2 + A' 2 /2) 2 - (A 4 - A' 4 )) 1/2 . 


(186 - /) 







Subsituting Eqs. (18) together with (14) into Eq. (8), the right-hand side of the wave 
equation becomes 


— / —(7 CI T iJ cy )e ,wt =‘brit/3 0 T ^-8(x)8(ij)e l “ z '' : " 

c J 0 2ir/u |e| 

[(K ± Tk q /2)a ± a ± e u!< ^ k T^ ,19) 

+(K*Zk q /2)t3?:a T e t(K ** k '' 2)z } . 


Dispersion R e la tion 

To obtain the dispersion relation the coupled differential equation must be solved, 
subject to the boundary conditions on the waveguide. To this end, we assume that the 
electron beam propagates within a perfectly conducting cylindrical waveguide of radius 
r g . In general, the complex amplitude for the right- and left-hand circularly polarized TE 
waves are written as, 

A ± (r,0,z)e~' ut = 6 ±nm [ < / n _ l (^ nm r)e ± + / n +i(/inm’')« ±,2i, eil 

exp [i (k± nm z ± (n - 1)0 - ud)], (20) 

where J n is the nth order Bessel function, b± nm are complex constants. are real 

constants determined by the boundary conditions and n = 1,2,3,... and m = 1,2.3.... are 
waveguide mode indicies. 

The boundary condition ia such that the tangential component of the electric field 
vanishes on the waveguide surface, r = r g . At r = r g we have E • e# = 0, where E is 
the total electric field and eg = — texp(— z’0)e + -f iexp( — i0)e_ is the unit vector in the 
azimuthal direction. Applying the boundary condition at r = r ? , we find the condition 
J' n {fAnm r g) = 0) so that n nm r g equals the mth positive zero of J’ n . 

Since A is driven by an effective transverse line current, we expect TE lm modes will 
be excited. Substituting the TE\ m mode representation into the left-hand side of Eq. (4). 
we obtain two sets of coupled equations. 
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E - Hi™ - «™) -UfW)4 ±lm 

E- Hi™ - *™) * 0 w!« ;w '^ 

4tr |e| / 2 "” dt | 

— o/ 0 / e (“T l *^y) 

c m 0 c 2 y 0 2717 u; 

Operating on Eq. (21) by both fq* d0 f^ 3 J a ((j. lm r)T<ir and 
/ 0 2ir dffe*' 20 Jg’ •MMim»*)rdr and solving for b± l<m , we find that 


"I.'*!, 

[■M/*imr)(J« T t^yjl^o + / 2 (Mimr)(7 t i tJy)l r=0 e*' 2 * 1 . 

( 22 ) 

where / lm = /'* [J£(Mi m r) + J?(/ti m r)] rdr = “ 1 WmA 2 ). 

The dispersion relation greatly simplifies when only one waveguide mode, say the 
TEn, takes part in the interaction, i.e., n = 1, m = 1, a ± = b ±u . k± = k~ n . 
Substituting (19) into the right-hand-side of (22), we obtain 

-fcs(fc±-^) I> ± 6 ±n e‘ fc ^ 

L V C - v o J ( J 3 V 

where fcjj = 2u/^ 0 Iu — 2v^i\ l l{^ 0 (^\ l r 2 - l)Ji(/*nr 9 )). The RH and LH circularly po¬ 
larized waves are simultaneously excited only when k+ = fc_ — k q . Eliminating 6 in . the 
dispersion relation coupling the RH and LH circularly polarized waves becomes 


[„ k;(k - ijjlv 0 +■ k q ) 2 D- 

1L 

,) 2 D*1 

[ u > 2 /c 2 - (k + k q ) 2 - n 2 u 

Jr -/< 2 -* 2 

- MIi 1 

- k*K* ( 7 

1 ( k - m / v a + k q 

L_) 

h ’ w/f-v-riJy 

vu i 2 ic 2 - (k 4 - k q ) : 

- Mil / 






where k = k +, A:_ = fed- A:,, A' = K± — k - u/v„ 4- k q / 2, £>± = A' 2 T A'i A' - A'f and 
R = R± = D+D _ — A” 4 = A 4 — (2A| 4- A’ 2 )A’ 2 4- A* 4 — A* 4 . Equation (24) can be put into 
the form 


D + 


kli* - u/v 0 + k q ) 2 i r _ ki(k-u/v 0 ) 2 
» 2 /c* — (As 4- k q ) 2 — Mii J [ ■ * 2 /c 2 - fc* - Mh . 


- A' 4 = 0, 


and/or 


A = 0. 


(25a) 


(256) 


Equation (25a) agrees with the dispersion relation in Ref. 9 with the vertical field set 
to zero. Equation(25b) is the dispersion relation of the hybrid cyclotron and quadrupole 
modes with image fields in a waveguide in the absence of electromagnetic fields. 

The dispersion relation possesses the following symmetry: 


{k,k q ,B a ) s=> (k - k q , -k q ,-B 0 ). 

Utilizing this symmetry condition, the discussions in Section III can also be applied to 
rotating quadrupoles of the opposite helicity. The polarization of the waveguide mode 
associated withe the wave number k that the three-wave is unstable is determined by the 
polarization of the waveguide mode that intersects the unstable beam modes. 

The modes taking part in the interaction can be conveniently classified by setting the 
beam current equal to zero. The dispersion relation, for zero beam current, reduces to 

[k 2 - a , 2 /c 2 + Mu] [(fc 4- k q ) 2 - <*i 2 j c 2 + Mu] 

[(* - u/vo + fc,/2) 2 - (d 2 + d 2 )] [(* - c/ Uo 4- k q /2) 2 - (d 2 - d|)] = 0, (26) 

where 

4 = ~(A„(A' 0 — k q ) 4- k q /2), 

4 = \(Kl(k 0 - k q f + 4K-k 2 ) 1 ' 2 . 

The first and second bracketed terms on the left-hand side of Eq. ( 26 ) represents the 
LH and RH circularly polarized transverse electric waveguide modes, while the third and 
fourth terms, R , are hybrid cyclotron and quadrupole modes. 
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A diagram of the dispersion relation in the zero beam current limit is shown in Fig. 
2. The parameters for the plot are guide magnetic field B a = — 1 kG, k q — 0.5 cm -1 , 
quadrupole gradient B q k q = 200 G/cm, j a = 5 and drift-tube radius r g = 3 cm. The 
curves are associated with the RH and LH polarized waveguide modes respectively. The 
cut-off frequency of the waveguide TEu mode is /i n c. The two straight solid (—) lines 
correspond to the beam lines associated with the d\ — d\ expressions, and the two dashed 
(- - -) lines correspond to the beam lines associated with the d\ 4- d\ expression. 
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III. STABILITY REGIMES AND ANALYTICAL EXPRESSIONS FOR GROWTH 
RATES 

The dispersion relation, Eq. (25a), contains i) a region of orbital instability (in the 
absence of the electromagnetic waves), ii) regions of three-wave instability, and iii) regions 
of stablility. In this section, we will obtain the conditions delineating the various regimes 
and find analytical expressions for the maximum growth rate in each of the three-wave 
unstable regions. The stability diagrams are obtained in (Je q ,K a ) space for given values of 
7 o, r g and B q k q in the limit of zero beam current. 

Orbit Unstable Regime 

The expression R = 0 is the dispersion relation for the particle dynamics in the pres¬ 
ence of stellarator windings with an axial magnetic field. This expression is in agreement 
with Eq. ( 10 ) of Ref. 2 in the limit of perfectly conducting walls. The electron beam in 
this configuration can be unstable when (d\ — d\) <0. The unstable values of K 0 are 

K e rit,2 = ^ - 2 K q <K 0 < A' e rtt,3 = | + 2 A'„ (27) 

where K q = fl 9 /v 0 , fi, = \e\B q /j 0 m 0 c, K a = fi 0 /v„ and Sl 0 = |e|B 0 / 7 0 m 0 c. Equation (27) 
is in agreement with the stability condition of Ref. 1 in the limit of straight cylindrici’ 
geometry and zero beam current. It is interesting to note that this condition is also in 
agreement with the condition for beam envelope stability in the limit of zero space charge . 13 

Three-Wave Unstable Regimes 

The three-wave instability will occur when, for example, the RHCP waveguide mode 
intersects, in the ( u>,k ) plane, the appropriate beam mode given by 

(k-v/v 0 + k q /2) 2 ~(d\ -dl) = 0, (2S) 

and {d\ - d\) > 0. For k q > 0, the instability for the RHCP waveguide mode occurs for 
<ji > 0. For K 0 < K e rit,2i three-wave is unstable (Region I) when the RHCP waveguide 
mode intersects the beam line w/v„ = (k 4- k q /2) + \f d\ - d\. For K 0 > A' cr , ( , 3 , the three- 
wave is unstable (Region II) when the RHCP waveguide mode intersects the beam line 
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u/v a = {k -f- kq/2) — yjd\ — d|. Identical three-wave instability growth rates occur for the 
LHCP waveguide mode, for u>/c with sign opposite to the unstable inodes associated with 
the RHCP waveguide mode. 

Three-Wave Stable Regime for K„ < A'-,,, •> 

Stability is achieved when the waveguide cut off frequency jz u c is sufficiently large so 
that intersection with either of the beam lines, defined by Eq. (28) cannot be achieved. 
The condition in terms of the waveguide mode cutoff is 

<ZMn > *<i + 2 (dl - dj) 1/2 , (29) 

where 



Based on Eq. (29), the region of 

fc, > qfin and K a < K C vit,z (30) 


is alway three-wave unstable. If the inequality in (29) can be satisfied, we can solve for 
the explicit value of K 0 for the three-wave stable regime. 

For k q < qnn and K 0 < K erxt< 2 , condition for stability in terms of the variable K 0 is 


2fK 0 (K 0 - M + f 2 - 4A*, k q < 0, 


where 

/ = qnu{k q - qnn/ 2 ). 

Defining 

C = fc*(l+8A, 2 /)-2/. 

we solve for K 0 with k q < qn n , and find four situations: 
i) for / > 0 and C, > 0, the stable range of K a is given by 


K c rit,l 



c 1/2 

— < Ko < K e rit.2' 


(31) 


(32) 


(34o) 
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ii) for / < 0 and C > 0, the stable values of K a are 


K„ < A' cr it,i = smaller 


°f[ 



iii) for / > 0 and ( < 0, all values of 


(346) 


K 0 < K cr 


(34 c) 


are unstable, 

iv) for / < 0 and ( < 0, all values of 


(34d) 


are stable. 

Three-Wave Stable Regime for K n > K,^ i 

The three-wave interaction is also stable when the RHCP waveguide mode intersects 
only the top (the beam line with larger u> for the same k) but not the bottom beam mode 
assocated with Eq. (28), and K a > K cr it, 3 - This occurs when 

qmi>k q -2{d\-d\f'\ (35) 


The three-wave interaction is stable for 

fc, < qnu and. A„ > A' cri t, 3 - (36) 

For kq > qn\\ and K 0 > A' er ,t ,3 , the values of K 0 that are three-wave stable are 

k C 1/2 

A'o > K C rit,4 = 4-(37) 

In this regime, one can show that f > 0 and £ > 0. In the limit of small quadrupole 
gradient, luge 7 <, and K a > A cr tt, 3 , the stability condition is approximately 11 


K 0 > fc, -mu/7o. 


(38) 
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The various operating regimes are illustrated as functions of k q and K a in Fig. 3, for 
7 „ = 5, r g = 3 cm and quadrupole gradient B q k q = 200 G/cm. This plot assumes that 
the quadrupole gradient is a constant, and {k q , K 0 ) are allowed to vary. The horizontal 
separation of the boundaries for the orbital unstable region is 4A' ? . Since the stability 
boundaries are obtained in the limit of zero beam current, the area of the actual stable 
regions will shrink slightly as the value of the current is increased. 

Analytical Expressions for the Growth Rates 

In the three-wave unstable regimes, we can obtain analytical expressions for the peak 
growth rates. The dispersion relation (25a) for the coupled RH and LH polarized wave 
can be rewritten as 

(a, 2 - u/ 2 )((a/ -a/ 2 ) 2 - Aa/ 2 )((u/ - a;,) 2 - Au/ 2 )(u/ 2 - u, 2 ) = 9, (39) 


where 


9 = k\v 2 0 c 2 [(u> - v„k) 2 (u 2 - wl)D+ + (a/ - v 0 (k + fc,)) 2 (u> 2 - w 2 )£>_] , 

t*>x = y/k 2 + Mii c i W 2 = v 0 (k 4- k q /2), a/ 3 = y/(k + k q ) 2 + jtjjc, Ao/J = u*(d? + d;), and 
Ao/| = t/ 2 (d? — d\). Based on the numerical results of the full dispersion relation in Eq. 
(25a), instability occurs at the intersection of the RHCP waveguide mode and one of the 
two modes of Eq. (28) for fc, > 0 and o/ > 0. Defining u = u>i 4- and <r = the 

dispersion relation reduces to 

(u, 2 _ u/ 2 )((u; - o/ 2 ) 2 - Au/ 2 )((o/ - o/ 2 ) 2 - Au/|) = (7, (40) 

where <r = k 2 h v 2 0 c 2 (w, - t».*) s 0+|, w±w 0+= (Au/ 2 /c) 2 ± A'i( Awj/c) - A', 2 
and the top and bottom signs in D+ refer to the three-wave unstable regions I and II 
respectively. Equation (40) is the dispersion relation when the RHCP waveguide mode is 
uncoupled to the LHCP waveguide mode. 

The instability region I in Fig. 3 is the result of the waveguide mode intersecting the 
upper mode (the beam line with larger w for the same k) given in Eq. (28). We will write 
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w = wj + Su and wi = u> 2 4- Aw 2 . We will assume Ao/i - Au> 2 >> 6w. The dispersion 
relation is given approximately by 


/Su Aw 2 \/£uA 2 _ <r/c 3 1 

\ C + C A C J " 2«1 (Aw 2 -Aw;)' 


(41) 


For simplicity, we will assume the temporal growth rate is much smaller than the separation 
of the beam modes, i.e., 6u << 2Aw 2 . With this assumption, we obtain a simple expression 
for the dispersion relation, 


/6w\ 2 ^ <r/c 2 1 

\ c / 4a?!Au) 2 (Awj — Awf) 


(42) 


For values of K 0 in the unstable region I, i.e., K 0 < K cr it, 2 , the values of cr are positive, 
and Eq. (42) gives the temporal growth rate. 

When the RHCP waveguide mode intersects both modes given by Eq. (28) with 
k q > 0, w > 0 and K 0 > K C rit,3 t the instability occurs only at the intersection of the 
RHCP wave mode and the lower beam mode. This gives the instability regime II in Fig. 
3. We will define W! = w 2 — Aw 2 . The dispersion relation becomes 


^ ^Au 2 j y _ 


1 


*/c 3 _ 

2wj (Aw 2 — Awf)* 


(43) 


Here again, we assume 2Aw 2 >> 6w and the dispersion relation reduces to 

/ Su \ 2 ^ <r/c 2 1 

\ c / 4 u/iAu/ 2 (Aw 2 — Awf) 


(44) 


For K C rit, 3 < K a < k q , the quantity cr is negative, and Eq. (44) gives the temporal growth 
rate in region II of Fig. 3. The analytical expressions (42) and (44) show that the temporal 
growth rate T = /m(£w) scales approximately as the square root of the beam current, i.e., 
r/c oc l}' 2 . 
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rv. NUMERICAL RESULTS 


The full dispersion relation, Eq. (25a) is solved numerically to: i) obtain the growth 
rates and group velocities, ii) verify the various operating regimes and the analytical expres¬ 
sions for the temporal growth rates and iii) show the scaling of the growth rate with respect 
to the various parameters. The numerical studies center around the parameters shown in 
Table I. With the quadrupole wave number k q chosen to be 0.5 cm' 1 (A, = 12.57 cm), 
we may demonstrate each of the different operating regimes by varying the axial magnetic 
field, 2?<>, except in one of the stable regimes, where we take k q =0.1 cm" 1 . 

Figure 3 is a plot of the various operating regimes in the parameter space of k q as a 
function of K 0 for 7 <> = 5, B q k q = 200 G/cm and r g = 3 cm as in Table I. 

Numerical Results from Full Dispersion Relation 

The dispersion diagram with current /(, = 1 kA for the five different regimes are shown 
in Fig. 4: 

a) B a — —1.0 kG ( K 0 = -0.24 cm' 1 ) and k q = 0.5 cm' 1 in the three-wave unstable 
region I, 

b) B„ = 2.15 kG ( K a = 0.26 cm' 1 ) and k q = 0.5 cm' 1 in the orbit unstable regime, 

c) B 0 = 3.5 kG (K 0 = 0.42 cm' 1 ) and k q = 0.5 cm' 1 in the three-wave unstable region 

II, 

d) B„ = 5.0 kG ( K a = 0.60 cm' 1 ) and k q = 0.5 cm -1 in the three-wave stable regime, 

e) B 0 = —5.0 kG ( K a = —0.60 cm' 1 ) and k q = 0.1 cm' 1 in the three-wave stable 
regime. 

Plots of the temporal growth rate as a function of wave number k are given in Figs. 5-7 
for each of the different regimes with current h = l kA. Figure 5 plots the temporal growth 
rates in region I, for B a = —1.0 kG ( K 0 = -0.12 cm' 1 ), B 0 = 0, B a = 1.0 kG ( K„ = 
0.12 cm -1 ), and B 0 = 1.3 kG ( K a = 0.156 cm -1 ). Only the growth rates associated 
with Rt(v) > 0 are plotted; these are associated with RHCP waves. The growth rates 
associated with Re(u>) < 0 are identical and are associated with LHC’P waves. The growth 
rate and the range of unstable values of k increase as K a approaches the orbit unstable 
value of A’erit.z = 0.154 cm -1 . 
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Figure 6 shows temporal growth rates in the three-wave unstable region II, for B 0 = 
2.95 kG (K a = 0.35 cm' 1 ), B 0 = 3.0 kG (K a = 0.36 cm' 1 ), B a = 3.25 kG (K„ = 
0.39 cm -1 ) and B 0 = 3.5 kG ( K 0 = 0.41 cm~ l ). For some values of B 0 , the RHC'P wave 
intersects the beam line two times as indicated by temporal growth rates at two separate 
regions of wave number k for a given axial magnetic field. The growth rate and the range of 
unstable values of k increases as K a approaches A' C Ht, 3 - The stability boundaries predict 
that the three-wave instability is present for K C rit ,3 = 0.35 cm -1 < K 0 < A' cr ,t ,4 = 
0.41 cm -1 (2.95 kG < B a < 3.5 kG). Numerical results for /*, = 1 kA are in excellent 
agreement with the theory. 

Figure 7 plots the temporal growth rate versus k in the orbit unstable region for 
B 0 = 2.15 kG ( K 0 — 0.26 cm' 1 ). The region of instability covers both positive and 
negative values of wave number k, even at regions where the beam lines are far from the 
waveguide mode. Thus, the unstable growth rate is the result of unstable beam orbits in 
the stellarator and guide fields, irrespective of the electromagnetic waves. 

As the beam current increases, the coupling of the RHCP and LHCP waves becomes 
stronger and the actual stability boundaries deviate from the analytic expressions. Figures 
8-10 are plots of temporal growth rate as a function of wave number k for I b = 10 kA. 
Figure 8 plots the temporal growth rates in region I, for B 0 = -1.0 kG (K a = —0.12 cm -1 ). 
B 0 = 0, Bo = 1-0 kG (K 0 = 0.12 cm” 1 ), and B a = 1.3 kG ( K a = 0.156 cm" 1 ). 

Figure 9 plots temporal growth rates versus k in the three-wave unstable region II, 
for Bo = 3.0 kG {K a = 0.36 cm -1 ), B a = 3.25 kG ( K a = 0.39 cm~ l ), B a = 3.5 kG (K a = 
0.41 cm~ 1 ) and B a = 3.85 kG ( K 0 = 0.45 cm -1 ). The range of B 0 for the three-wave 
unstable region II at zero beam current is 2.95 kG < B 0 < 3.5 kG. For current of 
I p = 10 kG, the RHCP and LHCP waves become coupled for B 0 = 2.95 - 3.25 kG. 
At the same time, the three-wave unstable regime is extended beyond B a - 3.5 kG to 
B 0 = 3.85 kG. The RHCP and LHCP waves become coupled for B a = 2.95 - 3.25 kG and 
the three-wave unstable regime is extended into a region where K 0 > Kent a = 0.41 cm _1 . 
where the analytic {I\, =0) expressions predict stability. 

Figure 10 plots the temporal growth rate versus k in the orbit unstable region for 
B 0 = 2.15 kG ( K 0 = 0.26 cm -1 ). Here the region of instability covers essentially all values 
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of wave number k. 

A summary of peak temporal growth rates as a function of normalized guide field K 0 
for h = 10 kA is shown in Fig. 11 for 7o = 5, 7 * = 7.5 and = 10. The group of curves 
on the left belong to unstable region I and the group of curves on the right belong to 
unstable region II. The gap separating the two groups of curves corresponds to the orbit, 
unstable region. Plots of the group velocity as a function of K a are shown in Fig. 12. 
Those values plotted are associated with the growth rates of Fig. 11. The group velocity 
approaches the beam velocity in three-wave unstable region II. 

We also examine the effect of the waveguide radius r g on the three-wave instability. 
As the guide radius is increased, the waveguide cut-off frequency fi\\c decreases. For the 
three-wave unstable region I, the intersection of the waveguide mode and one of the beam 
lines in Eq. (28) can occur for negative values of wave number k. When this occurs, the 
pha»e velocity is negative and group velocity is positive but reduced in value. Figure 13 
plots the temporal growth rate and the group velocity as a function of K 0 for r g = 5 cm 
and = 10 kA, while keeping ail the other parameters the same as in Table I. 

CpmpfrrigQfl at Analytic and Numerical, Rggnlts 

The analytical expressions for the temporal growth rates for the two regions of the 
three-wave instability are given in Eq. (42) and (44). Figure 14 is a comparison of the 
maximum temporal growth rate as a function of the normalized guide field A' 0 for the 
values obtained from the numerically solved frill dispersion relation (solid curves) and 
from the analytical expressions (dashed curves) for current h = 1 kA. The agreement 
between the analytical and numerical results of the temporal growth rates is good as 
longhair Awj/c > F/c (see Eqs. (39)-(44)). On the dispersion diagram, 2 Au> 2 /c is the 
difference in £requency/c of the beam lines (28) that could go unstable. To illustrate the 
case Awj/c > T/c, we take K„ = —0.36 cm -1 (B a = —3 kG ). The numerical result of the 
temporal growth rate is T/c = 4.3 x 10 -3 cm -1 , which is much smaller than Au ; 2 /c = 0.2 r >. 
The analytical result for the temporal growth rate is also T/c = 4.3 x 10 ~ 3 cm' 1 . As Au): ; c 
decreases and T/c increases, the analytical expressions for the growth rate becomes less 
accurate. At K 0 — U.12 cm' 1 ( B 0 = 1 kG), Au> 2 jc = 8.2 x 10 -2 , which is comparable to 
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T/c. The analytical expression for the growth rate, T/c = 2.50 x 10 2 cm 1 , is 8% larger 
than the numerical result of T/c = 2.30 x 10 -2 cm -1 . 
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V. CONCLUSIONS 


The addition of stellarator windings to an axial guide field was proposed as a method 
of transporting a high current beam in a curved geometry with a high tolerance to energy 
mismatch. 1 The stabilty properties of such configurations have been clarified in the present 
study, which included the beam centroid motion, the electromagnetic waveguide inodes, 
expressed in terms of right-hand circularly polarized (RHCP)' and left-hand circularly 
polarized (LHCP) waves, and the induced image forces on the conducting boundaries, in 
addition to the external magnetic fields. 

We find five operating regimes: i) two physically distinct three-wave unstable regimes, 
ii) an orbit unstable regime, in which the beam centroid is unstable independent of the 
electromagnetic waveguide modes and iii) two stable regimes. We have obtained analytical 
expressions for the boundaries of the various stability regimes in parameter space and have 
presented algebraic expressions for the growth rates in each of the two three-wave unstable 
regimes. These analytical results are valid in the limit of low beam current, when the 
RHCP and LHCP waves decouple, and are in good agreement with those obtained via 
numerical solutions of the full dispersion relation. 

The simplified dispersion relation for uncoupled RHCP and LHCP waves, given in 
Eq. (40), has been shown to produce results that are in close agreement with those of the 
full dispersion relation. The simplified dispersion relation is not valid, however, when the 
coupling between the RHCP and LHCP waves is strong, as in the orbit unstable regime. In 
three-wave unstable region II, the coupling is strong only at high current. The growth rates 
in this regime, for example, from the full (coupled) dispersion relation are about 20 — 30% 
larger than from the uncoupled dispersion relation for k w = 0.5 cm -1 and / = 10 kA. In 
three-wave unstable region I the RHCP and LHCP waves are essentially uncoupled. 

These results suggest that the three-wave instability can be avoided by appropriately 
choosing the various parameters. Results show that as the beam energy increases, the 
stability conditions become more restrictive. In such cases, it may be necessary to decrease 
the quadrupole gradient, B q k q , quadrupole wave number k q and/or increase the value of 
the guide field B 0 in order to remain in the stable regime. 
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Table I: Parameters Used iu Section IV 


Quadrupole gradient, B q k q 
Beam energy (j 0 ) 

Beam current, 

Drift-tube radius, r g 

Calculated parameters 

TE n cut-off frequency/c, 
qn 11 

Kqk q 

& = vjc 

For rotating quadrupole wavelength. A 7 = 12.57 cm 

Wave number, k q 
K q 

K C rit,2 {Bo = 1.32 kG) 

Kcrit'3 {Bo = 2.95 kG) 

Kent,* {Bo = 3.52 kG) 


200 G/cm 
5 

1,10 kA 
3 cm 

0.614 cm ~ 1 
0.256 cm' 1 
0.024 cm -2 
0.9793 


0.5 cm -1 
0.048 cm" 1 
0.154 cm -1 
0.346 cm -1 
0.413 cm -1 
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g the coiulurting cylindrical drift tube, the rotating quadn 








Wave Number, k (cm ') 

Fig. 2 The dispersion diagram in the limit of zero beam current. The curves are the LH 
and RH circularly polarized TEn waveguide modes. The dashed (- - -) and solid (—) 
straight lines are the beam modes. 









t of (k Vi K„) parameter space showing the various operating regimes. Stability 
indaries are calculated from Eqs. (27)-(37) for parameters of Table I. 
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Fig. 5 Plots of the temporal growth rate versus k in three-wave unstable region I. for 
h = 1 kA and 

a) B 0 = -1.0 kG ( K 0 = -0.12 cm -1 ), 

b) Bo = 0, 

c) Bo = 1.0 kG ( K 0 = 0.12 cm" 1 ) and 

d) B 0 = 1.3 kG ( K 0 = 0.156 cm' 1 ). 
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Fig. 6 Plots of the temporal growth rate versus k in the three-wave unstable regie 
current = 1 kA and 

a) B a = 2.95 kG {K„ = 0.35 cm -1 ), 

b) B a = 3.0 kG (A'„ = 0.36 cm' 1 ), 

c) B a = 3.25 fcG (A„ = 0.39 cm -1 ) and 

d) B a = 3.5 kG (A' 0 = 0.41 cm' 1 ), 
















Wave Number, k (cm 1 ) 

Fig. 8 Plots of the temporal growth rate versus k in three-wave unstable region I, foi 
/(, = 10 kA and 

a) B 0 = -1.0 kG (K„ = -0.12 cm' 1 ), 

b) B a = 0, 

c) B a = 1.0 kG ( K 0 = 0.12 cm -1 ) and 

d) B 0 = 1.3 kG (K 0 = 0.1.56 cm' 1 ). 











Temporal Growth Rate/c, T/c 



Wave Number k (cm" 1 ) 

Fig. 9 Plots of the temporal growth rate versus k in the three-wave unstable region II. fo 
current I\> = 10 kA and 

a) B a = 3.0 kG (K 0 = 0.36 cm- 1 ), 

b) B 0 = 3.25 kG ( K 0 = 0.39 cm" 1 ), 

c) B a = 3.5 kG (A'o = 0.41 cm- 1 ) and 

d) B 0 = 3.85 kG (A'„ = 0.45 cm' 1 ). 
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Fig. 10 A plot of the temporal growth rate versus k in the orbit unstable region, for current 
I h = 10 kA and B a = 2.15 kG ( K a = 0.26 cm' 1 ). 
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in Table I. 
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Fig. 14 Comparison of the maximum temporal growth rate as a function of A' 0 for the values 
obtained from the numerically solved full dispersion relation (solid curves) and from 
the analytical expressions (dashed curves) for the parameters in Table I with current 
I h = 1 kA. 
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